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ASPECT  RATIO  AND  LOADING  EFFECTS  ON  A  CRACKED  RECTANGULAR 
TENSION  SPECIMEN  MADE  FROM  A  LINEAR  ELASTIC  MATERIAL 

D.T  Baron 
Raytheon  STX  Corp. 

Edwards  AFB,  C A 

C.T.  Liu 

Air  Force  Research  Laboratory 
Edwards  AFB,  CA 


ABSTRACT 

Finite  element  calculations  are  performed  to  determine  the  effect  of  specimen  aspect  ratio  (specimen  height  {h) 
divided  by  specimen  width  (w) )  on  the  specimen  geometry  stress  intensity  correction  factor  {F  =  JjE/{oJt^)) 
(where  J  is  the  J  -integral,  E  is  Young’s  modulus,  a  is  the  remote  force  divided  by  the  uncracked  area,  and  a  is  the 
crack  length).  F  is  a  function  of  specimen  geometry,  and  modifies  the  value  of  the  Mode  1  stress  intensity  (F^) . 

Linear  elastic  fracture  mechanics  (LEFM)  and  plane  stress  are  assumed.  Two  specimen  types  are  simulated,  center 
cracked  tension  (CCT),  and  single  edge  notched  tension  (SENT).  Two  boundary  conditions  are  simulated,  prescribed 
uniform  load,  and  prescribed  uniform  displacement.  For  each  of  the  four  combinations  of  specimen  type  and 
boundary  condition,  an  estimate  is  made  of  the  minimum  value  of  h/w  for  which  the  effect  on  F  can  be  neglected, 
in  which  case  F  can  be  calculated  using  an  equation  instead  of  from  a  numerical  analysis. 

INTRODUCTION 

An  important  problem  in  engineering  design  is  evaluating  the  strength  and  reliability  of  structures.  It  is  known 
that  the  strength  of  a  structure  may  decrease  during  the  period  of  its  design  life  under  service  loads.  One  of  the 
common  causes  of  weakening  is  the  development  of  structural  cracks.  However,  the  presence  of  cracks  in  a  structure 
does  not  necessarily  mean  that  the  structure  is  no  longer  functional.  Therefore,  to  calculate  the  ultimate  strength  or  the 
maximum  service  life  of  a  structure,  studies  are  conducted  to  quantify  the  significance  of  flaw  types  and  size,  and  to 
determine  whether  a  crack  will  grow,  and  the  rate  of  growth.  In  this  study,  finite  element  simulations  are  performed  of 
two  commonly  used  fracture  mechanics  test  specimens,  CCT  and  SENT  (Figure  1),  It  is  important  to  be  able  to 
calculate  the  Mode  1  stress  intensity  factors  for  these  specimens  for  any  combinations  of  boundary  condition, 
specimen  aspect  ratio  {h/w)  and  crack  length  to  specimen  width  ratio  {a/w) .  The  specimen  geometry  stress 
intensity  correction  factor  (F)  enables  the  calculation.  For  large  values  of  h/w ,  F  is  independent  of  h/w ,  and 
formulas  exist  for  F .  In  this  paper,  finite  element  calculations  are  made  to  determine  how  big  h/w  must  be  in  order 
to  be  neglected,  and  the  effect  on  F  for  smaller  values  of  h/w . 

SPECIMEN  GEOMETRY  STRESS  INTENSITY  CORRECTION  FACTOR 


The  specimen  geometry  stress  intensity  correction  factor,  (F) ,  is  defined  as 


F(y,  F,  a,  a)  = 


(1) 


where  J  is  the  /-integral,  E  is  Young’s  modulus,  a  is  the  remote  force  divided  by  the  uncracked  area,  and  a  is  the 
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crack  length.  F  quantifies  the  effect  of  specimen  geometry  on  Ky.  Equation  (1)  has  the  form  of  a  calculated  value  of 

Ky  (JjE) ,  divided  by  a  reference  value  of  Ky  (a^/7Ca) .  If  for  a  particular  experimental  specimen,  F ,  a ,  and  a  are 
known,  then 


(2) 

This  is  how  F  is  used,  for  calculating  the  values  of  Ky  for  experimental  specimens.  For  any  LEFM  specimen  with  a 
remotely  applied  load,  Ky  can  be  written  in  the  form  [1] 


(3) 


P  is  an  applied  concentrated  force,  b  is  the  specimen  thickness,  and  w  is  the  specimen  width.  Each  specimen  type 
has  a  particular  function  f{a/w) .  For  CCT  and  SENT  specimens,  the  validity  of  equation  (3)  depends  on  the 
specimen  aspect  ratio  (h/w)  being  large  enough  so  that  their  exists  a  cross  section  between  P  and  the  crack  which 
is  loaded  in  uniform  tension.  In  that  case  the  particular  value  of  h/w  has  no  effect  on  the  value  of  Ky,  Substituting 

the  identity  P/{bJw)  =  ajw  into  equation  (3)  gives 


(4) 


Equation  (4)  can  be  used  with  the  LEFM  relation  J  =  Ky/E  to  transform  equation  (1)  into 


(5) 


Note  that  equation  (5)  depends  only  on  the  ratio  a/w .  Geometrically  scaling  a  specimen  (multiplying  a,  h,  and  w 
by  the  same  scalar)  has  no  effect  on  F,  Using  equation  (5)  with  the  function  f{a/w)  for  a  CCT  specimen  [1]  gives 
the  CCT  specimen  correction  factor, 


^CCT  - 


(6) 


Using  equation  (5)  with  the  function  f(a/w)  for  a  SENT  specimen  [1]  gives  the  SENT  specimen  correction  factor, 


(7) 


Equations  (6)  and  (7)  are  the  CCT  and  SENT  specimen  geometry  stress  intensity  correction  factors  for  large 
specimen  aspect  ratios  (h/w) ,  The  values  of  ^sekt  computed  for  any  valid  value  of  a/w 

(0  <  a/w  <  1)  using  equations  (6)  and  (7).  No  numerical  analysis  is  required. 
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NUMERICAL  RESULTS 

Table  1  shows  the  values  of  finite  element  calculations  of  the  specimen  geometry  stress  intensity  correction  factor 
(F) .  The  calculations  were  made  to  determine  the  influence  of  specimen  aspect  ratio  {h/w)  on  F  for  both  CCT  and 
SENT  specimens.  In  the  table,  each  simulated  specimen  is  labeled  as  being  of  type  CCT  or  SENT.  Additionally  each 
specimen  has  a  three  character  code  consisting  of  an  uppercase  letter  (“L”  or  “D”)  and  a  two  digit  number.  “L” 
indicates  the  specimen  was  subjected  to  a  uniform  load,  and  “D”  indcates  the  specimen  was  subjected  to  a  uniform 
displacement.  The  two  digit  number  is  the  specimen  aspect  ratio.  For  instance,  specimen  “CCT  LI 4”  is  a  simulated 
CCT  specimen  with  a  prescribed  uniform  load,  and  h/w  -  1/4 .  There  are  two  exception  specimens  in  Table  1.  The 
correction  factors  shown  for  specimen  “CCT  EQ6”  were  calculated  using  equation  (6),  and  those  for  specimen 
“SENT  EQ7”  were  calculated  using  equation  (7).  No  finite  element  calculations  were  made  for  these  two  specimens. 
The  values  of  F  in  the  row  for  specimen  “CCT  EQ6”  are  the  limiting  values  for  a  CCT  specimen  as  h/w  increases, 
and  the  values  of  F  in  the  row  for  specimen  “SENT  EQ7”  have  the  analogous  meaning  for  SENT  specimens.  Figures 
2-5  were  made  from  the  data  of  Table  1,  and  are  plots  of  (F)  versus  a/w  for  constant  h/w .  Figures  2  and  3  are  for 
CCT  specimens,  and  Figures  4  and  5  are  for  SENT  specimens.  Figures  2  and  4  are  for  applied  uniform  load,  and 
Figures  3  and  5  are  for  applied  uniform  displacement.  The  data  of  Table  1  specimens  “CCT  L##”  and  specimen  “CCT 
EQ6”  are  plotted  in  Figure  2.  The  data  of  Table  1  specimens  “CCT  D##”  and  specimen  “CCT  EQ6”  are  plotted  in 
Figure  3.  The  data  of  Table  1  specimens  “SENT  L##”  and  specimen  “SENT  EQ7”  are  plotted  in  Figure  4.  The  data  of 
Table  1  specimens  “SENT  D##”  and  specimen  “SENT  EQ7”  are  plotted  in  Figure  5.  Figure  7  plots  the  same  data  as 
Figure  5,  but  a/w  and  h/w  are  switched.  Figure  6  shows  second  order  polynomial  curve  fits  of  the  induced  nodal 
stresses  along  the  tops  of  specimens  “SENT  D##”  (Figures  5,  7).  Curve  fit  values  are  shown  instead  of  the  actual 
values  because  a  finite  element  solution  gives  actual  nodal  stress  values  which  oscillate. 

Figures  2  and  3  show  that  for  a  CCT  specimen  subjected  to  either  prescribed  uniform  load  or  prescribed  uniform 
displacement,  the  effect  of  h/w  on  F  can  be  neglected  when  h/w>  3  ,  and  equation  (6)  can  be  used  to  calculate  F . 
Figure  4  shows  that  for  a  SENT  specimen  subjected  to  a  prescribed  uniform  load,  the  effect  of  h/w  on  F  can  be 
ignored  when  h/w  >  2 ,  and  equation  (7)  can  be  used  to  calculate  F ,  Figure  5  shows  that  for  a  SENT  specimen 
subjected  to  a  prescribed  uniform  displacement,  it  is  always  necessary  to  include  the  effect  of  h/w  on  F  (except 
when  a/w  is  small  (say  a/w  <  0.10 ),  and  simultaneously  h/w  is  not  small  (say  h/w  >  3 ));  and  therefore  it  is 
necessary  to  perform  a  numerical  analysis  in  order  to  obtain  the  value  of  F  for  a  particular  value  of  h/w . 

Figures  2-5  show  the  same  trends  for  both  specimen  types  (CCT,  SENT).  For  a  specimen  with  a  prescribed 
uniform  load,  decreasing  the  specimen  aspect  ratio  (h/w)  increases  the  correction  factor.  For  a  specimen  with  a 
prescribed  uniform  displacement,  decreasing  the  specimen  aspect  ratio  decreases  the  correction  factor.  Explanations 

can  be  given  for  these  trends.  The  equation  for  the  correction  factor  is  F  =  JJe/{<5J^)  .  For  the  case  of  prescribed 
uniform  load  (a  is  constant)  with  a  crack  of  fixed  length  (a  is  constant),  only  J  can  vary.  As  the  specimen  aspect 
ratio  decreases,  the  amount  of  material  between  the  (constant)  applied  traction  and  the  crack  surface  decreases.  The 
remaining  material  must  deform  more  to  carry  the  same  load.  This  results  in  a  larger  crack  opening  displacement 
(COD),  therefore  a  larger  value  for  ,  and  consequently  a  larger  value  for  J ,  resulting  in  a  larger  value  for  F .  For 

the  case  of  prescribed  uniform  displacement,  as  the  specimen  aspect  ratio  decreases,  the  amount  of  material  between 
the  (constant)  applied  displacement  and  the  crack  surface  decreases.  The  resistance  of  this  material  to  displacement 
decreases  (the  material  acts  like  a  beam  with  a  shrinking  moment  of  inertia)  while  the  resistance  of  the  material  in 
front  of  the  crack  is  constant.  The  induced  stress  distribution  at  the  top  of  the  specimen  rotates  so  that  the  stress 
behind  the  crack  tip  decreases  with  respect  to  the  stress  in  front  of  the  crack  tip.  This  results  in  a  decrease  in  the 

quantity  Jj/ a ,  and  therefore  a  decrease  in  F .  Figure  6  shows  the  approximate  induced  stress  distributions  along  the 
tops  of  the  SENT  specimens  subjected  to  a  prescribed  uniform  displacement  (specimens  “SENT  D##”)-  As  the 
specimen  aspect  ratio  (h/w)  decreases,  the  stress  distribution  rotates  counterclockwise  and  the  crack  is  unloaded 
(the  crack  is  on  the  left  side  of  the  specimen).  For  a  CCT  specimen,  any  cross-sectional  stress  distribution  must  be 
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symmetric.  This  makes  the  specimen  geometry  stress  intensity  correction  factor  (F)  for  a  CCT  specimen  subjected 
to  a  prescribed  uniform  displacement  much  less  sensitive  than  the  corresponding  SENT  specimen,  to  the  effect  of 
specimen  aspect  ratio  (h/w) . 

In  Figure  5  (SENT  specimens,  prescribed  uniform  displacement),  by  increasing  h/w  above  the  maximum 
shown  value  (8/1) ,  it  should  be  possible  to  produce  an  F  curve  which  is  arbitrarily  close  to  the  (limit)  curve  of 
equation  (7).  In  Figure  7  (a  re-plot  of  data  from  Figure  5),  the  limiting  value  for  F  as  h/w  increases  of  (for  instance) 
the  curve  for  a/w  =  0.5  ,  should  be  2.827  (the  value  of  equation  (7)  when  a/w  -  0.5 ).  By  inspection  of  Figures  5 
and  7  it  seems  that  h/w  would  have  to  be  very  large  for  the  effect  of  h/w  on  F  (especially  for  large  values  of  a/w ) 
to  be  negligible. 


DETAILS  OF  FINITE  ELEMENT  CALCULATIONS 

The  finite  element  calculations  were  made  using  the  ABAQUS  commercial  finite  element  program.  Every 
number  in  Table  1  (except  for  the  rows  corresponding  to  “CCT  EQ6”  and  “SENT  EQ7”)  required  one  finite  element 
calculation  (one  ABAQUS  run).  So  there  were  a  total  of  203  runs.  For  each  run,  the  mesh  and/or  boundary  conditions 
had  to  be  modified.  Instead  of  using  a  commercial  graphical  mesh  generator,  a  Fortran  90  program  was  written  which 
reads  a  file  containing  the  values  of  problem  definition  parameters  and  mesh  definition  parameters,  and  then  generates 
an  ABAQUS  input  file.  This  results  in  a  consistent  and  precise  method  of  mesh  generation.  For  each  specimen  type 
(CCT  or  SENT),  a  particular  mesh  (node  coordinates  and  element  -  node  connections)  was  generated  from  the  values 
of  six  parameters;  the  specimen  width  (w ),  the  specimen  height  {h ),  the  crack  length  to  specimen  width  ratio 
(a/w) ,  the  number  of  element  sectors  (n^) ,  the  number  of  element  orbits  (n^) ,  and  the  ratio  of  the  specimen  width 
to  the  radius  of  the  inner  element  orbit  (w/r^ ) .  Figure  8  shows  examples  of  the  meshes  that  were  used.  The  program 
attempts  to  make  all  of  the  element  sector  angles  equal.  The  curve  of  the  inner  element  orbit  is  always  circular. 
Assume  for  each  mesh  that  the  origin  of  a  coordinate  system  is  located  at  the  crack  tip,  the  x  axis  is  parallel  to  the 
crack  with  positive  direction  pointing  toward  the  specimen’s  uncracked  edge,  and  at  the  y  axis  is  perpendicular  to  the 
crack  with  positive  direction  pointing  towards  the  top  of  the  mesh.  Along  each  of  the  three  axis  sector  lines  (+;r  ,-x, 
+y )  the  growth  rate  of  the  distances  between  adjacent  element  orbit  lines  (the  orbit  space  growth  rate)  is  constant. 
The  growth  rate  along  one  of  these  sector  lines  is  calculated  from  the  sector  line’s  length,  the  number  of  element 
orbits,  and  the  radius  of  the  inner  element  orbit  (rj ) .  For  all  203  runs,  four  of  the  six  mesh  parameters  were  constant, 

i.e.,  w  =  1 ,  =  48  ,  Wq  =  48 ,  w/r^  =  1000 .  h  was  varied  according  to  the  specimen  aspect  ratio  {h/w) ,  and 

{a/w)  was  also  varied.  This  means  that  the  same  basic  mesh  was  used  for  all  203  runs.  ABAQUS  input  files  could  be 
quickly  generated,  and  in  one  instance  77  were  made  within  an  hour.  Eight  node  rectangular  isoparametric  plane 
stress  elements  were  used  with  full  integration.  Each  of  the  crack  tip  elements  was  made  triangular  by  locating  three 
nodes  along  one  element  side  at  the  crack  tip.  All  of  the  crack  tip  nodes  were  tied  together.  The  side  nodes  of  the 
crack  tip  elements  were  moved  to  the  quarter  points.  48  calculations  of  the  /  -integral  were  made  for  each  ABAQUS 
run  (each  element  orbit  line  was  used  for  a  contour).  For  almost  every  run,  the  J  values  varied  for  only  a  few  of  the 
innermost  and  outermost  contours.  Usually  the  J  values  for  the  bulk  of  a  particular  run’s  contours  were  identical  to 
each  other.  It  so  happened  that  for  the  runs  when  a/w  =  0.01 ,  the  orbit  space  growth  rate  along  the  -x  sector  line 
(the  crack  sector  line)  was  less  than  one.  This  is  because  >  a/48,  i.e.,  w/1000  >  0.01  w/48.  (This  means  that 
was  bigger  than  the  average  of  the  48  orbit  spaces  along  the  crack  sector  line  (if  equaled  a/48  ,  then  the  orbit 

space  growth  rate  would  equal  one).)  However,  for  these  cases  Table  1  seems  to  show  that  the  calculated  J  values  are 
still  reasonable.  Each  SENT  specimen  mesh  models  one  half  of  the  symmetric  specimen  (dividing  the  specimen 
along  the  crack).  Therefore,  the  aspect  ratio  for  a  SENT  specimen  mesh  is  equal  to  h/{2w) ,  i.e.,  it  is  one  half  of  the 
aspect  ratio  for  the  corresponding  SENT  specimen  {h/w) .  In  Figure  8,  all  of  the  meshes  are  for  SENT  specimens,  so 
the  aspect  ratio  of  each  of  the  three  shown  meshes  is  equal  to  one  half  of  the  aspect  ratio  of  the  associated  specimen. 


In  Figure  8,  for  each  mesh  the  location  of  the  crack  tip  is  visible  at  the  bottom,  and  the  crack  is  along  the  bottom  to  the 
left  of  the  crack  tip.  Each  CCT  mesh  modeled  one  quarter  of  the  bi-symmetric  specimen  (dividing  the  specimen  along 
the  crack  and  across  the  crack’s  midpoint).  So  the  aspect  ratio  for  a  CCT  specimen  mesh  is  equal  to  the  aspect  ratio 
for  the  same  CCT  specimen.  For  both  types  of  specimen  meshes  (CCT,  SENT)  the  mesh  bottom  to  the  left  of  the 
crack  tip  (the  crack)  was  subjected  to  a  zero  traction,  and  so  was  the  right  edge  of  the  mesh.  The  mesh  bottom  to  the 
right  of  the  crack  tip  was  constrained  from  movement  in  the  y  direction.  The  top  of  the  mesh  was  constrained  from 
movement  in  the  x  direction.  The  top  of  the  mesh  was  either  subjected  to  a  uniform  displacement  or  a  uniform  load  in 
the  y  direction.  For  SENT  specimen  meshes,  the  left  edge  was  subjected  to  a  zero  traction.  For  CCT  specimen 
meshes,  the  left  edge  was  constrained  from  movement  in  the  x  direction.  Table  1  indicates  that  25  of  the  203  runs 
failed.  All  of  the  failed  runs  were  for  a/w  <  0.10 .  In  these  cases  ABAQUS  complained  that  some  element  areas  were 
negative.  It  is  unlikely  that  any  element  area  was  literally  negative;  i.e.,  it  is  unlikely  that  one  part  of  an  element’s 
boundary  crossed  another  part.  “Negative”  probably  means  a  mathematically  prohibited  shape  (For  a  simpler 
element,  a  four  node  rectangle,  an  example  of  a  prohibited  shape  is  when  one  of  the  interior  angles  is  greater  than  n ). 

CONCLUSIONS 

The  specimen  geometry  stress  intensity  correction  factor  F  -  4Je/{(5J^)  is  independent  of  specimen  size 
(multiplying  a,  h,  and  w  by  the  same  scalar  doesn’t  change  ^ ).  ^  is  dependent  on  specimen  proportion 
{a/w,  h/w) .  For  both  types  of  specimens  considered  (CCT,  SENT),  for  the  boundary  condition  of  prescribed 
uniform  load,  decreasing  the  specimen  aspect  ratio  {h/w)  increases  F .  For  both  types  of  specimens,  for  the 
boundary  condition  of  prescribed  uniform  displacement,  decreasing  h/w  decreases  ^ .  For  any  typical  SENT 
specimen  subjected  to  prescribed  uniform  displacement,  the  effect  of  h/w  on  F  cannot  be  neglected,  and  for  each 
set  of  geometry  parameters  {a/w,  h/w)  a  finite  element  calculation  is  required  to  determine  F .  For  SENT 
specimens  subjected  to  prescribed  uniform  load,  the  effect  of  h/w  on  F  can  be  neglected,  and  equation  (7)  can  be 
used  to  calculate  F ,  when  h/w  >  2 .  For  CCT  specimens  subjected  to  either  prescribed  uniform  displacement  or 
prescribed  uniform  load,  the  effect  of  h/w  on  F  can  be  neglected,  and  equation  (6)  can  be  used  to  calculate  F ,  when 
h/w  >  3 . 
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SPEaMEN* * ** 

a  /  w  =  0.01 

a  /  w  =  0.05 

a/w  =  0.10 

a  /  w  =  0.20 

a  /  w  =  0.30 

a  /  w  =  0.40 

a  /  w  =  0.50 

SENT  L14 


SENT  L21 


SENT  L31 

** 

SENTL41 

1.147 

SENT  L51 


SENT  LSI 


SENT  D14 


SENT  D13 


SENT  D12 


SENT  Dll 

0.950 

SENTD21 

1.110 

SENTD31 

*s)e 

SENT  D41 

1.037 

SENT  D71 


SENT  D81 


Table  1.  Finite  element  values  of  the  specimen  geometry 
stress  intensity  correction  factor,  — . 


*  Each  specimen  is  indicated  by  a  type,  either  “CCT”  or  “SENT”,  an  upper  case  letter,  either  “L”  or  “D”,  and  a  two 
digit  number.  “L”  indicates  that  uniform  load  is  prescribed.  “D”  indicates  that  uniform  displacement  is 
prescribed.  The  number  indicates  the  specimen  aspect  ratio.  For  example,  specimen  “SENT  D21”  is  a  SENT 
specimen  with  prescribed  uniform  displacement,  and  h/w  =  2/1 . 

**  No  value  determined.  ABAQUS  error  message  said  that  some  elements  indicated  negative  areas  for  calculations 
at  some  integration  points. 


CCT  EQ6  ==>  Values  for  CCT  specimen  are  calculated  using  equation  6.  (No  finite  element  calculations.) 
SENT  EQ7  ==>  Values  for  SENT  specimen  are  calculated  using  equation  7.  (No  finite  element  calculations.) 
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Figure  4.  Finite  element  values  of  correction 
factor  vs.  a/w  for  constant  h/w. 

SENT  specimens  with  prescribed  uniform  load. 


equation  7 
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hfw  =  8/1 


.A/f =.1/1 
hfw  =  1/2 


Figure  5.  Finite  element  values  of  correction 
factor  vs.  a/w  for  constant  h/w. 

SENT  specimens  with  prescribed  uniform  displacement. 


